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ABSTRACT 

Simulation as a mode of computer use in instruction 
has been neglected by educators. This paper briefly explores the 
circumstances in which simulations are useful and presents several 
examples of simulation programs currently being used in high-school 
biology, chemistry, physics, and social studies classes. One program, 
STERIL, which simulates screw-worm fly population changes with the 
use of pesticides and sterilization techniques, is discussed in 
detail. <RB) 
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DIGITAL SIMULATION IN EDUCATION 



I. Introduction 

Educators have found many ways of using computers. At the high-school 
level , computers are used primarily for administrative purposes while at the 
college level, the computer is used for educational purposes as well as for 
administrative ones. These educational uses at the college level have been 
associated primarily with programming courses, or with large-scale or diffi- 
cult computations (e • g • , in nonlinear, or distributed systems) , 

One mode of computation which has been largely ignored among educators 
is the simulation mode. One notable exception is Huggins 1 JOBSHOP (1). Med- 
ical schools are beginning also to explore the uses of simulations in devel- 
opment of diagnostic technique (2), as well as in other ways. 

Hugging (3) has characterized the digital computer as a highly- flexible 
laboratory in which we may explore any environment we wish. It is with this 
mode of computer utilization that we shall he concerned in this paper. 

The author has been working for several years with high— school teachers 
in explorations of the utility of computer simulations in their curricula. 
Simulation programs have been found to be effective in biology, chemistry, 
earth science, physics, and the social sciences (4) . 

In the paragraphs below, we shall explore the circumstances under which 
simulations are useful, and shall present several examples of simulation 
programs which have been developed recently. 
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XI , A Rationale for Simulation 



Few people would dispute the desirability of learning by direct experience. 
The investment of billions of dollars by our universities, and by our primary 



and secondary schools in laboratory facilities is evidence enough. Generation 
after generation of students re-measure the acceleration due to gravity, or 
re-discover Ohm’s Law, or study the giant chromosomes of the drosophila fly, 
even though these things have been done by very capable scientists employing 
the most elaborate facilities available. How much easier (and how much less 



expensive) it would be to describe Millikan’s oil-drop experiment, and its re- 
sults, along with pictures of the apparatus, in a textbook — yet many physics 
laboratories in high schools and colleges are equipped with this apparatus, and, 
every year, legions of students repeat this experiment. Their e/m ratios are 
very imprecise compared to the carefully-determined value of the physicist; how- 
ever, they have benefited from having had the experience themselves. This xs as 
true for non-scientists, who never again will have a laboratory experience, as it 
is for the physicist who will spend his life in the laboratory. 



Even though direct experience is very desirable, our students frequently 
are prevented from such learning experiences. Some of the reasons for this are 

1) The necessary equipment is not available because of expense 

or it is too complex or delicate to permit students to 
use it (e.g., in high-energy physics) . 

2) The sample size available in the real,. world is too small to . . 

permit generalizations (this .is especially , true in the 

training of medical students in the diagnosis of disease, 
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Medical students in Hew York State, for example, run 
into very few cases of malaria, and develop little 
experience with it) s 

3) The experimental technique is difficult and must be developed 
over an extended period (e.g*, in experiments in genetics, 
and in titration) . 



by There are serious dangers to the student (e.g. , where radia- 
tion or high temperatures are involved, where there may 
be explosive mixtures of gases, or where highly-toxic 
materials are required) * 

5) The time scale is too short or too long to permit tne student 
to make observations (e.g., the runaway of a nuclear 
reactor — here, of course, there are other reasons 
for not permitting students to do the experiment — or 
the study of the dynamics of populations) * 



6) The opportunity to experiment directly is not available 

(e.g* , in studies of economic, politeal, or social sys- 
tems, or in studies of human genetics, or spread of 
disease)* 

7) When it is desirable to measure variables which are diffi- 

cult to access (e.g., the tension on a pendulum string, 
or the differential effects of the gravitational effects 
of the earth and the moon on an orbit ting satellite) . 

8) When measurement and other noise obscures the important 

phenomena (in; the computer, we can create a world in 
which there is ho ; noise and in which instruments are 
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perfect ^ and then show the student how these things obscure 
the data of interest) s 

9) There are times when it is useful to underscore the significance 
of natural laws (e . g „ , study of non-inverse-square-law 
gravitational systems (5) , or non-Hendelian genetics). 

If any of the foregoing circumstances exists, the student normally is 
excluded from learning by direct experience. Typically, here, he learns by 
reading a textbook, or by attending a lecture on the subject. In both cases, 
the learning Is less than' completely satisfactory. (Who would accept a sur- 
geon who has learned surgery solely from books and lectures?) 

Simulation is an educational tool which can overcome the foregoing in- 
hibitions; and which can provide, in many cases, a far better learning situa- 
tion than can be provided by either textbooks or lectures. 

It should be emphasized that simulation is proposed as a complement to 
existing educational techniques, rather than as a substitute for any of them. 

III. Examples of Educational Simulations 

Several examples of simulations developed by the author and his colleagues 
are described below. These simulations were developed for high-school students 
in biology, chemistry, physics, and social studies; however, they are appropriate 



for college students, at least at the freshman level. They have been used 
successfully by such students at two: colleges (the Polytechnic Institute of 
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Aj_STEB.il 

One pest-control technique which has proven to be very success- 
ful in control of screw-worm flies is the "sterile-male" technique, 
in which male flies are raised in a laboratory, sterilized by irradia- 
tion, and, then, released in infested areas. These sterilized males 
compete with normal males for the females. The success of the tech- 
nique depends upon the fact that females mate only once. Those which 
mate with sterile males lay unfertilized’ eggs, and fail to reproduce. 
This technique is described in an interesting paper by Knipling (6), 
and in an extensive report (7) . 

A simulation program called STERIL (developed by A. Irishman , 

C. Losik, and this author) permits the student to explore the relative 
merits of pesticides, sterilized males, and combinations of these two, 
in the control of screw -worm flies. The student is able to develop, 
and try out, various policies for controlling, and eliminating, these 
pests. 

Several sample runs of this program are shown in Figs. 1-7, In 
Fig. 1, a pesticide is applied every seven days, but no sterile males 
are released . Here, there is a steady— state oscillation with a seven- 
day period and an average level of 0.25 million normal males. The 
population increases after each application of pesticide, because of 
immigration into the area, and because of eggs laid by earlier gen- 
erations, - '■ • ' v -y r.-'-j. , . : : 

In Fig. 2, pesticide is applied every day for 24 days only (again. 



me 



no sterile males are released) , On the twen th- f i f tli day, the popula™ 
tioxi begins to grow again* because of immigration; and, by the seven* 
tie th day , it re turns 
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The approach shown in Fig. 3 consists of the application of pesti- 
cide only once on the first day. From then on, one-hundred-thousand 
sterile males are released each day. Here, the normal population de- 
clines continuously beginning on the sixth day; until, on the seventy- 
fifth day, it is reduced to sixty thousand. 

Fig. 4 shows the result of a policy similar to that of Fig. 3, 
except that only ten-thousand steriles are released each day. This 
policy decreases the steady— state normal population only to 880,000. 

One interesting characteristic of this policy is the dips at days 
21, 41, and 61, The twenty-day period is a result of the twenty-day 
period from the laying of eggs to the maturation of the new females 
and their laying of eggs . 

The result of application of pesticide once, and release of 100,000 
sterile ma les each day for twenty-four days is shown in Fig. 5. Again, 
because of immigration and maturation of already— laid eggs, the popu- 
lation returns to 1,000,000 sterile males. Here, the steady-state 
average is 0,18 million compared to 0.25 million in Fig. 1, This 
slight reduction, in spite of the large number of sterile males re- 
leased each day, occurs because the pesticide is non-sel ective. Sterile 
and normal males are affected equally. 

The policy of Fig, 7 (pesticide only for the first ten days, and 
then no control) shows that early termination of control causes more- 
rapid return to normal population levels than with ex tended control 
programs -- either pes ticide or sterile males (c.f . , Figs . 2 and 6) • % r 
Clearly, this program . permits exploration of a circumstance which 
the student could not study directly, except perhaps in graduate school. 



and which he could not study analytically because of the complexity 
of the equations, 

B. KINET 



The program KINET permits the student to explore the dynamics 



of a chemical system in which a specific compound is converted to 



a single daughter- compound , The student specifies the forward-rate 
constant, and the equilibrium constant. From these parameter value 
the computer produces plots of the concentrations of mother and 
daughter species as functions of time. 

Four sample runs of this program are shown in Figs, 8-11, 

KINET was developed by J, Marchisotto and H. Shannon, 



C. SLITS 

One of the classic experiments in the elucidation of the wave 
properties of light was Young's Double-Slit Experiment* This ex- 
periment can be performed directly, even by high-school students; 
however, because of equipment limitations, there are severe re- 
strictions on parameter ranges, and on the depth of exploration 
available to the student, SLITS is a simulation of this experi- 
ment which permits the student to explore parameter values over 
a wide range, to observe, in detail, the variation of light in- 
tensity along the viewing screen (this is extremely difficult 
to instrument in the laboratory) , and to draw inferences about 
the inter-relationships among the variables-. /' 

Three sample runs of this program are shown in Fig. 12. 

This program was developed by A. Caggiano. 



-a- 



D. GENE1 

In the study of genetics, students usually have difficulty appre- 
ciating the probabilistic nature of the transfer of genetic traits 
from parents to offspring* GENE1 has as its purpose the development 
of an understanding of the fact that the Mendelian ratios apply only 
when large numbers of offspring are involved. 

Sample runs of this program are shown in Figs. 13-16 . In these 
runs, there are increasing numbers of offspring. Theoretically, for 
this set of genetic traits, the Mendelian ratios are: genotype ratio 
— homozygous dominant: heterozygous: homozygous recessive - 1:2:1; 
phenotype ratio — dominant; recessive = 3:1. (In this example, 
brown eye color Is dominant over blue* 

For 10 offspring (Fig* 13), the actual ratios are far from the 
theoretical ones (1:6:3 and 2*33*1; rather than 1:2:1 and 3:1). As 
the number of offspring Increases, the ratios approach the theoretical 
values. When there are 10,000 offspring, the actual ratios are 
1:1.93:1.03 and 2,84:1 — within 5% of the theoretical values. 

This program was developed by the author. 

e * scatter ' ' /h 

This program simulates the phenomenon called Rutherford, or Coulomb 
scattering. Here, a series of energetic particles is fired at a tar- 
get particle . The output of one r (in of this :prograni is shown in 

Fig. 17. The target nucleus is located at the T in the center of , ;• yi 

the plot. There are -five equally-energetic particles released at 

different locations with, respect to the target particle. These £ 
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five particles scatter by an amount which varies inversely with 
distance from the target axis* * 

This program was developed by A. Caggiano and D * Scarl . 

F. EC ONI 



This program is a very simple-minded simulation of the economy 
of the United States* The plot of Fig. 18 shows the variations 
in net national income with various taxing, governmental spending, 
and investment policies, (This program was developed by the author.) 

IV* Conclusion 

Several examples of educationally-oriented simulation programs have been 
described. All of them have been found to be useful pedagogically ; and at least 



two (SLITS and STERIL) have been received very enthusiastically by students* 
One fact which has emerged clearly from this developmental effort is the 
need for a diversity of backgrounds and capabilities in the development of 
these simulation programs . 

The discussion of STERIL above is extensive primarily to show the versa- 



tility and relatively-high fidelity of the simulation model used here. This 



'.V . 

? •• * • 



program required the combined skills of a very capable. programmer, a knowledge- 
able bibldgis t , and a system engineer wi th experience in model building . , 






m 









' 'V-* fT-'. -V. • /. T 1 ■ 



iv. .V.- 



It is clear to this author that , before simulation begins to have the 

^ ,, - 

dramatic impact on education which it offers potentially, more subject-matter 

• ; ‘ ^ : \ V . : , • : • ' ‘ :■ :• . • ; -v V; - 

specialists .system engineers , and computer scxentxsts from the univarsities , 

^ s'/v r : v : ;• - •■•'S' r..; ; v -e ' ' - r ^ V v l V; -v : v.-- ~ ■ : ; • ~ • \ : 

. •- - - . - . • s • ‘ '• _ -- •— 

l from industry, must begin contributing to the development of material. 

. .. . . 



and 
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LET F « THE FORWARD HATE CONSTANT . 

USTK« THE EOUlLUHilUN CONSTANT FOR THE REACTION A=P 
TVI’E IN THE CONSTANTS F AND K IN iHAT ORDER 
? 5/ .5 ' . 






PERCENT CONCENTRATION OFAC*) AND HC+) 
0 25 50 ... 75 100 • 

I--- I —'—I _I 

TIME - ’ * . ' 

Q + * 

.046 * 

.002 ' * * 

,13b + * 

. ItS A + • .• * 

.23 ■ + . * 

.276 + * • 

.322 , + ■ * 

* 36b * : +. * 

,414 + * 

.46 . + . * 

Fig. 8. KINET run with forward rate constant - 5; equilibrium constant =0.1 

, LET F - THE FORWARD liATB CONSTANT. . . :'i. 

LETK= THE EO JI LI Bril UN CONSTANT FOn T n E t H E A C T I 0 N . A = P 
TYPE IN THE CONSTANTS F AND K IN THAT ORDER - ; 

7 5 ,. 2 . . \ - " • 

; ’’ v;; 4Tf- •• * * * * *' * v.* ^ « *•* * * *'* ♦ * f « * * * * * * *' . • 

£a'..V ■■ /•. ■; ■ ' -> 

CONCENTRATlONliOFAC *>~>and p < + ) 
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LET F " THE FOWLED RATE CONSTANT 

LET K = THE. E2UILI PR I UN CONSTANT FOR THE REACTION A -P 
TYPE IN THE CONSTANTS F AND K IN THAT ORDER. 

? 5 * 0.1 



***! M =**'*********** ****** ****** 



PERCENT CONCENTRATION OF ACT-) AND PC + > 



TIME 

0 

.01254545 
.02509091 
•03763636 
• 0501 81 82 
.06272727 
*07527273 
.08781 fi IS 
. 100 3636 
.1129091 
• 125454 5 
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100 

-I 

* 
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Fig. 1C. 



KINET run with case forward rate ron«t»nt as in Fig > 
equilibrium constant - u.o 
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LET F 
LET K 



THE FONARD RATE CONSTANT A • 

LET K - THE EOUILIBHIUN , CONSTANT FOR THE REACTION 
type IN THE CONSTANTS F AND K INf THAT ORDEH. ,;TF:,;: : :: .^ 

WiT;?? Y’:= • „ tr, n "i- 

• V - v - . . : . rv .„ . . +#+**+************.**+*******+*■•■-_ - ; . . .. . . . 
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TOUNG'S DOUBLE SLIT EAPEitiMENf 
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ViiAf I b iTf:£ MEW SLI T SEr’AiuVl' iOM CD) IN M 1 LLI ME THiiS ? U*2 
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D - • 2 M 1 LL 1 M E 
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- • 26 » ♦ 

- ,24 * * 

- . 22 • * 

- , a • * 

- , ib . * 

-.16 . * 

- .14 • * 

- . 12 • * 
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WOULD YOU LI KB TO T*tY AMO 1‘nEn VALUE 
fchA i' IS THE NEW SLIT SEPARATION CD) 



OF D Cl-YEB.. U - Pi 0 3 ? 
IN I ¥ i I LL 1 METEitS ? I .0 
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Fig .12. continued . 

■ ■ V. ‘I- ■ • .■r.Y.'r 'l ^ - n : i h 

•.'*• -.-v. ’’ -•• : :V ; . • : 

'—•r : -'••• ■x •' ■" - T 

' • • . . • • ■ ' 

V- r ! ‘. c > 



M /Y/uY v: Y r: 




c. Same as Fig. 12. a. 
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WrlA r tUlE Ti:i.E TWO T.iAlXB TO OK STUDIED? 
DOL< 1 MAN r Tit A I T ? U;iO V;N’ 
nEUEiBlVE TAAIT? BLUE 



*****!* . 

\ 

GSNOTYBE OF FEMALE BAuENT? BiiOWtNLBLUE 
GEN OT V BE OF MALE B ABENT? BBOWM/ BLUE 
HOw MAN* OFFBBHING DO YOU WANT TO 3IUDX? 10 
DETAILED nEBQrfTC YE3 Oii NO)? TEB 



OFFSBiilNG NO 
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GEMOTiBE BAT TO r 6 s -v ;3' V: : T: ^Tyr'V " 

BnENOTTBE BATIO ' . : 3.333333*5 1 T*;:;'* /; /if- v-*;; ; 
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GENOTYPE OF' FEMALE PABENT? BuOWN* BLUE 



GENOTYPE OF KALE PA BE NT? BBOWN*6LUE 

how many offsbuing do you want to study? loo 

DETAILED. BEPOHTC YES OB NO'? NO 



2jij(e .+ K- + * * *■ V * =♦= * 

GENOTYPE BAT 1 0 1 : 1 .56620 7 : . S6B0 6u 

PHENOTYPE BAT 10 3 S 1 

v Fig, 14. GENE! output with 100 offspring. 



GENOTYPE OF FEMALE PABENT? BBQ WNjBLUE 
GENOTYPE OF MALE PABENT? SHOWN* BLUE: 

HOW MANY OFF SP Til NG DO YOU WANT TO STUDY? lUUU 



RATIOS' ONLY WILL BE T'YPED> BECAUSE OF 
LABCS NQ. OF QFFSBBI \ 



GENOTYPE xiATiO • 
PHEMO T .PE jiAT 1 0 



1 : <2 .Vi 1666? : 1 .US 

2 . ySttPSA S 1 



i#tf #####> mV## 



Fig. 15, GENE1 with 1,000 offspring. 
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GENOTYPE OF MALE PABEN Y? ,DBO : .vN*BLuE 

• . • . ' , ■ ' ■ ■ 

un:.i :■/. a \\ v A m-itAUA 1 AIR f)0 YOU WANT TO STUDY? lUpUO 
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